Functional renormalization group approach to the dynamics of first-order phase 
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We apply the functional renormalization group theory to the dynamics of first-order phase tran- 
sitions and show that a potential with all odd-order terms can describe spinodal decomposition 
phenomena. We derive a momentum-dependent dynamic flow equation which is decoupled from the 
static flow equation. We find the expected instability fixed points; and their associated exponents 
agree remarkably with the existent theoretical and numerical results. The complex renormalization 
group flows are found and their properties are shown. Both the exponents and the complex flows 
show that the spinodal decomposition possesses singularity with consequent scaling and universality. 
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First-order phase transitions (FOPTs) are ubiquitous 
in our diverse physics world. Basically, their dynamics 
is described by two distinct mechanisms [l[: one is nu- 
cleation and growth and the other is spinodal decom- 
position. FOPTs show discontinuity at their transition 
points at which the first (and higher) order derivatives of 
the chemical potential of the two phases involved change 
discontinuously. This is in sharp contrast to continuous 
phase transitions that have singularity at their transi- 
tion points [2j and are generally described by the elegant 
renormalization-group (RG) theory [3|, Consequently, 
there is a challenging question of whether FOPTs can 
also show any singularity and consequent scaling and uni- 
versality. Recently, it has been found 0, HI that the spin- 
odal decomposition is possibly controlled by instability 
fixed points (IFPs) whose exponents equal those of the 
Yang-Lee edge singularity Q. Further evidence for the 
existence of the IFPs was studied in the two-dimensional 
Potts model Q- However, more studies on these new 
and imaginary fixed points are still needed in particular 
to search support from other theories and methods. 

A possible way to study the existence of the IFPs is the 
functional RG (FRG) theory, which, like the perturbativc 
RG theory also originates from Wilson's idea of RG [3] 
but provides us both perturbative and nonperturbative 
aspects. The FRG theory provides a unified picture to 
describe a lot of phenomena in many different fields Q. 
It has also been applied to understand nucleation in the 
scheme of Langer's theory In these studies for 
a final potential Uk f ~ rri^ (j> 2 /2 + 7fc/^ 3 /6 + A/ C/ 4 /8, 
which is identical with its initial form, where TOfe, is 
the mass and jk f and A^, are the couplings in the fi- 
nal scale kf, a ratio R(h) with h = 9Afc / m? hk IS m ~ 
troduced to determine the region of validity of the ho- 
mogeneous nucleation theory. When Afe./mfe, is fixed, 
the spinodal line is approached from h — > 0. This cor- 
responds to the region in which m/ c/ and Afc, tend to 
zero simultaneously and hence the potential is left only 
with the term jk f 4> 3 /6- Although the nucleation theory 



described by the FRG theory does not valid in this re- 
gion, it still indicates the spinodal decomposition has a 
potential of </> 3 like with a massless mode, which is just 
the viewpoint in [J]. In the FRG perspective, we thus 
expect that the potential of the spinodal decomposition 
contains all odd order terms of the field <f> for symme- 
try reason, in comparison with the critical phenomena 
which contain all even-order terms with Z2 symmetry. 
To confirm this argument, we need to find the expected 
IFPs in the FRG scheme, determine the corresponding 
exponents, and then compare them with extant results. 
In the FRG scheme, one needs to solve the exact static 



Wetterich equation 11] with appropriate approximations 
to determine a fixed point and its exponents. We shall use 
the Blaizot, Mendez-Galain, and Wschebor (BMW) ap- 
proximation 12j, which contains full momentum depen- 



dence and has yielded excellent results for static critical 
exponents among others [13]. As the spinodal decompo- 
sition is kinetic in origin, we have to consider dynamics. 
However, dynamic flow equations have only been derived 
with the derivative expansion method |l4j . So, we shall 
first derive in the BMW scheme a dynamic flow equation 
as well as the decoupled static BMW equation from a 
dynamic Wetterich equation. 

Consider the Hamiltonian in the bare scale A, 
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The dynamic action S for Model A [15( can be expressed 
in the form [la ]. 

S = I {^(d^/dt + Ddn^/S^) ~ D^ 2 } , (2) 

J X 

with a response field 4>, where x = (t,x) and D is a 
kinetic coefficient. To transform to the FRG scheme, we 
add to S a term 



Rk = DR k 



1 

1 



(3) 
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with the static regulator Rk(q) — Z k (k 2 — q 2 )Q(k 2 - q 2 ) 
[l7j . where k is the flow scale, $ = ((/>,<}>), Zk is the 
renormalization factor for the field <fi, the Heaviside 
step function, and r denotes a transposition. AS k plays 
a role of mass that suppresses fluctuations of momentum 
\q\ < k. As a result, the scale-dependent effective action, 



Tk[<P,<p] 
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with W fc [M] = In/ 2?$exp(-5-A<S fc + / x £' =1 Mi)> 
interpolates between the mean field action S[ip] = T k= \ 
and the effective action T[ip] — T k= o as fluctuations are 
progressively taken into account as k flows from the bare 
microscopic scale A to the long-distance scale 0, where 
ip = (tp, (p) with ipi denoting the average of 4>i and (h, h) 
is the conjugate source. We have introduced a subscript 
i 6 {1,2} such that (4>x,4>i) — (0,0) for convenience. 
Defining the 2x2 matrices of the vertex functions as [l4[ 



f^/V Y ^ - S2Tk 

L k V 1 ) 2 J T 7 7 r • 



\(n>2) 



(xi, • 



5 (n-2)f( 2) (x ls x a ) 



<Vi 3 (x 3 )---c>^ in (x„)' 
one writes the flow of the effective average action as 
f 



(5) 



(6) 



d s T k = -Tr / G k d s R k 



in the Fourier space, where G k = [T ( k > + R ki 
full field-dependent propagator, q = (ui,q), and s = kd k 
is the RG time. Equation (J7J is just the dynamic Wet- 
terich equation [8j] reflecting the RG flow of T k under 
an infinitesimal change of k. From Eq. (JT)) and using 
Eq. one can derive the flow equation of two-point 
vertex function in uniform field configurations 18 1 



is the 



0.rg?(p) = Tr / (d s R k )G k (q) 



1 



-lr$(p,-p,q) 



+rf ) (p,q)G fc (p + q)rf (-p,p + q)j G fc (q), (8) 

where are the matrix elements of . In Eqs. ([5]) 
to ([5]). all the ipi dependence are implicit. 

In order to derive from Eq. © the specific dynamic and 
static flow equations for Model A, one needs considering 
only those vertex functions that comply with its struc- 
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ture. Since only T 



.(3,1) 



= -DAa, and T 



(Li) = ^ 

(0,2) 



-2D exist at </?i = in 
the bare scale, we set other vertex functions zero [l9[ in 
Eq. ((SJ. To decouple the flow equations for dgT^ip) an d 

(2) 

the real part of d s T\ J(p), we drop the cj dependence in 
T k °' 2 ^ and assume 



,(i,i) 



(q) 



(0,2) 



(9) 



where T, n \q) is the static vertex function. Equation © 
is a kind of Ward identity expressing the fluctuation- 
dissipation theorem. It becomes exact when the uj de- 
pendence of F k °' 2 ^(g) is re-included and no mode cou- 
plings are present [20[. By Eq. ©, vertex functions of 

(2 1) f3 1) 

higher orders such as and 1^ ' depend no longer 

on lu cither. In addition, the external frequency is set 
to zero, since we mainly focus on p = (0,p). After 
performing the trace and integrating over uj, one finds 

(2) 

the static Wetterich equation about 1^ (p) decoupled 
from the dynamic behavior. Then, the approximation 

r *T • • • , gn, 0) = (?!,••• ,q n ) results in the 

BMW equation 



d s r k 2) (p)= I \dsRk(q)]G 2 (q) 

J a 



x \G k { P + q ) [a„if >(p)] - -d 2 r< 2) (p)| , (io) 

and a momentum-dependent dynamic flow equation 



d s vr>{ P ) = / [a s i? fc (g)]^( g )G fc ( P + g )rr ; (g) 

"'9 

xlf' 2 >(p + g) [9 v r( 2) (p)] 2 (A + 2 J B)(A + i?)- 2 , (11) 

(2) 

which is one of our main results, where Gk(q) = (9) + 



(7) R k (q)]-\ A ee Tf 2 \q)G- k 1 {p + q), and B = T k °> 2) (p 



q)G k ~ 1 (q). Note that the dynamics part does not con- 
tribute to Eq. (|10[) and the fixed points are thus solely 
determined by the latter as should be. 

Having these flow equations, one needs to drop their 
dimensions to access the fixed points. The dimensionless 
renormalized variables and functions are defined as 



p = p/k, q = q/k, 



U k = k d U kl 0=[k- 2 - d Z k ]V 2 ip, D = Z n D 



'k = 
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vf 2) = -pZoZkD}- 1 ^ 

(12) 

where U k is the potential, K^ 1 = d2 d - 1 Tr d / 2 Y{d/2) {Y 
is the Euler Gamma function), and Ze> and Z k are the 
renormalization factor of D and tp, respectively, and are 

1/2 — 1 /2 

related through Z k = Zr>Z k . In the scaling regime, 
Zd ~ k~ eh and Z k ~ k~ Vk , with e k and being related 
to the dynamic and the static critical exponents by z = 
2 — £ k= o and 77 = r?fc = o, respectively. Consequently, the 
dimensionless dynamic Sow equation satisfies 



<p=[k?- d Z k Kj 1 ]V*<p, 
<p, D = Z D 1 



(0,2) 



d s vf 2 \p) = ( Vk - e k +pd p ) ff 2 \p) + d s r k °' 2) (p), 

(13) 

where the last term is Eq. ([lip in its dimensionless form. 
The static flow equation can be separated into two di- 
mensionless flow equations: One is the flow equation of 
U k = ri ; (0) (a prime denotes a derivative with tp) ob- 



tained by setting p 
flow equation of \k 



in Eq. flTDJl, the other is the 
Z k (tp)/Z k (0) which comes from 



3 



d s z k = [d(d s rW(p))/d P * 



They are [21 1 



0.056 



d s U' k ' = (r, k - 2)U' k ' + -(d-2 + r, k )^U' k " 



1 - 



Vk 
d + 2 
1 



Giur-2Gl(Un 2 ,(14) 



1 - 



>7fe 
d + 2 



4G£E«-Gfttf (15) 



for our choice of i?fc, where G k = (1 

From the dimensionless dynamic and static flow equa- 
tions (p~3|) to (flSl) . one can find the IFPs and their ex- 
ponents in a way similar to the critical phenomena. The 
only difference is that we need to change to the symmetry 
that describes the spinodal decomposition. We use the 
field expansion method to find the IFPs 21 , 22j . The key 
point is how to retain terms in U k and Z k . As mentioned 
above, the potential should be of the form 



Uk = A 2 <^ 2 + A 3 ^ 3 + A 5 <^ 5 + • • • + \ ni $ ni + 



(16) 



Xk = 1 



■ X4<f 4 ■ 



Xn 2 <P 



(17) 



Eq. (fT3| recovers the dynamic flow equation in Ref. [14 1 
if the (p dependence of the latter equation is neglected, 
which confirms ours. From Eqs. (|15[) and (|13[) . one finds 



Vk 



Zk 



Glxl-±GlU' k " X 'k + {U> k ") 2 G 



Xk + dT2 G kXk - d+2 G WkXk 



Vk 



Vk 
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WVGk 



(18) 



(19) 



= 



n H 



d + 2-r] 
d+2--r) + 2z 



d-2 



d + 2~r) + 2z 



(20) 



which reflect the hysteresis of the coercivity and the rem- 
nant magnetization, respectively. 




o.ooo 



Re(^) 



The expansion of Z k (or Xk) is delicate. Any of its in- 
dependent term will bring even besides odd order terms 
to the unrenormalized potential. In this potential, how- 
ever, the even-order terms that arise from the even-order 
terms of Z k are generated by the coupling of the latter 
terms with A2 only, which vanishes at the mean-field or 
bare-scale spinodal point. Consequently, in order to have 
a potential of odd-order couplings at the bare scale, we 
retain only the even-order terms in Z k and write 



Substituting Eqs. (Hp and (17]) into Eqs. §14$) and (jIS]> . 
one gets a set of nonlinear flow equations for the coef- 
ficients X ni and Xn-2 ■, which reach their fixed points at 
d s X ni = d s Xn 2 — 0. To obtain the exponents, we use 
p = and (p — as the renormalization point and let 
ffe ' 2 ^?) = 1 in Eq. (fTT|) . With this approximation, 



respectively. In order to compare with direct numerical 
results, we also compute the hysteresis exponents [1,0], 
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FIG. 1. (Color online) (a) Stability of Im(A3) with increasing 
N in d = 3,4,5. RG flows of (b) Re(?7 fc ) and Im(A 3 ) and 

(c) Im(?7fc) and Re(As) from s = to s = —12 with four 
different values of Re(A°): 0.1 x 10" 4 , 0.8 x 10~ 5 , 0.6 x 10 -5 , 
0.4 x 10 -5 (top to bottom). Whereas the flows in (b) do not 
depend on Re(A-j), those in (c) are proportional to Re(Ajj) and 
can thus overlap when scaled by the proportionality, (d) The 
solid line and the dashed line are the maximum of Im(?7fc) and 
Re(As), respectively. Both of them are linear with Re(A|j). 
This relation is broken with larger Re(Ajj). In (b), (c), and 

(d) , d = 3 and TV = 9. 



We now present the results for the fixed points. For 
d > 6, the only infrared-stable fixed point found is a 
Gaussian fixed point with A ni = 0, confirming the mean- 
field behavior in these dimensions. For d < 6 on the 
other hand, a lot of new fixed points appear. To find 
the nontrivial IFPs, we start with n\ = 3 and «2 = 0. 
There are only three fixed points: Besides the Gaussian 
one, the others are just the IFPs of a pair of two purely 
imaginary conjugate values of A3. This can be checked 
by their stability against increasing n\ and n 2 . When 
ni is fixed, the IFPs become stable with increasing 712, 
whose maximum value is n\ — 3 for closing the set of 
equations. So, the stability can be checked by n\ — 2N + 
1 = 3, 5, • ■ • and 712 = ni — 3. One sees from Fig. QJa) 
that the imaginary part of A3, denoted by Im(A3), indeed 
stabilizes quickly as N increases. 

Having determined the IFPs, we can then extract their 
exponents. We note first that if we neglect the depen- 
dence of Xk on {p from the beginning, we find r\ = —0.176, 
—0.427, and —0.669 for d = 5, 4, and 3, respectively, 
which are already quite good compared with those ob- 
tained from a perturbative RG theory [23j shown in Ta- 
blelU Taking into account the field expansion of Xk spoils, 
however, the results. Yet, if we set X2 — in Eqs. (fT8]) 
and (fT9|) , the resultant exponents given in Table [I] agree 
remarkably with the extant ones. Two remarks are in 
order here. First, the perturbative results were derived 
only to two- and three-loop orders and thus may not be 
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TABLE I. Instability Exponents: a This work. b Perturbative 
RG Results [H, c Numerical results 



d 


6 


5 


4 


3 


v a 





-0.147 


-0.328 


-0.517 


v b 





-0.147 


-0.328 


-0.524 


z a 


2 


1.922 


1.809 


1.661 


z b 


2 


1.938 


1.874 


1.809 


n°k 


2/3 


0.650 


0.636 


0.624 


n c H 


0.655(2) 


0.645(3) 


0.625(12) 


0.595(30) 


" m 


1/3 


0.260 


0.168 


0.0546 


"m 


0.33(1) 


0.27(4) 


0.17(8) 





accurate. Second, our zs are a bit smaller than the oth- 
ers, this may arise either from our setting T^' 2 \q) = 1 in 
Eq. (fTTj) . or our omitting the (p dependence in Eq. (flUl) . 
Still, the remarkable agreement in all the three dimen- 
sions considered indicates that the condition of \2 = 
is not accident. A consequence is that it leads to a rela- 
tion between A2 and A3 via rjk from Eqs. (fl"4f and (fl"5l) . 
the relation which may be pertinent to the fact that the 
ip 3 theory has only one independent static exponent. Al- 
though further studies are clearly needed, the good agree- 
ments still confirm our identification of the imaginary 
fixed points with the IFPs. 

To show how to reach the imaginary IFPs from the 
bare scale, we fix the value of A2 to its fixed point value 
at any scale, and set Im(A° i>3 ) = 10~ 16 , Re(A" 1>3 ) = 0, 
and Xn 2 = at the bare scale indicated by the super- 
scripts, while varying Re(A 3 ). Then, we solve the set of 
nonlinear flow equations of A ni and Xn 2 U P to N = 9 
in d — 3 to 5. The flows start from the bare scale 
A (s = 0) to the zero scale (s — > —00), see Fig. QJb) 
and (c). The flows of Re(r]k) and Im(A 3 ) do not depend 
on Re(Ag); they overlap completely for different Re(A 3 )s. 
The flows of Im(?7fe) and Re(Aa), on the other hand, de- 
pend on Re(A 3 ). Nevertheless, they are proportional to 
Re (A3). In other words, if Re(A 3 ) is multiplied by a con- 
stant, then the flows can still overlap after multiplying all 
values of the flows by the same constant. If the value of 
Re(A 3 ) is too large, however, the proportional relation is 
broken [see Fig.[IJd)], and the IFP becomes unreachable. 
Although this appears to indicate that the IFP has a con- 
trolling region, it becomes larger as Im(A" i ) increases. In 
addition, flows of purely real A? s cannot reach the fixed 
points in our numerical analysis, though those from an- 
alytical solutions may. Therefore, we see that the IFPs 
are indeed infrared stable and reachable. A tiny initial 
complex coupling appears sufficient for the flows to reach 
the fixed point in the FRG theory. This confirms a pre- 
vious perturbative analysis [Ijj. As the complex flows 
connecting the physics at the bare scale with the physics 
at the zero scale are a process of adding fluctuations to 
the mean-field theory, which is controlled by the Gaus- 



sian fixed point and describes the spinodal decomposition 
well for d > 6 0, |j| , the FRG theory thus shows that the 
same phenomena in d < 6 should be controlled by the 
IFP, though the initial purely imaginary coupling may 
have to be identified. 

In summary, we have applied the FRG theory to the 
dynamics of FOPTs. The momentum-dependent dy- 
namic flow equation in the BMW scheme has been de- 
rived. We have found the IFPs with their exponents in 
good agreement with extant results albeit with some sub- 
tleties. This gives a physical meaning of spinodal de- 
composition for a potential with all the odd-order terms 
as compared with the critical phenomena for a poten- 
tial with all the even-order terms. The complex flows 
and their properties have also been shown in the FRG 
scheme. Both the exponents and the complex flows show 
that the spinodal decomposition has a behavior of singu- 
larity and consequent scaling and universality. 

This work was supported by the NSF of PRC 
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